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Abstract. A variational technique to describe the ground and scattering states 
below the break-up threshold for a three-nucleon system is developed. The method 
consists in expanding the wave function in terms of correlated Harmonic Hyper- 
spherical functions suitable to handle the large repulsion contained in the nuclear 
potential at short distances; three body forces have also been considered. The inclu¬ 
sion of the pp Coulomb repulsion in the p-d processes does not cause any particular 
problem, since no partial wave decomposition of the interaction is performed. Accu¬ 
rate numerical results are given for ground state properties and scattering lenghts, 
phase shifts and mixing parameters at three different energies of the incident nu¬ 
cleon. The agreement with other available results and with experimental analyses 
is higlly satisfactory. 

1. Introduction 

A sophisticated variational technique to describe the bound state of a three- 
nucleon system has been developed by the authors in ref.[l]. The wave function 
is expanded in channels, as in the Faddeev technique, and for each channel, the 
radial amplitude is, in turn, expanded in terms of correlated functions, constituting 
the Pair correlated Hyperspherical Harmonic (PHH) basis; such a basis results to 
be well suited to carefully take into account the correlations induced by the large 
repulsive terms of the nucleon-nucleon (NN) potential. This technique was applied 
in [1] to calculate the bound state w.f. of the triton, with the Argonne 14 (AV14) 
model of NN interaction [2], with results in complete agreement with those given 
by the best available techniques [3-5]. 

It is well known that focal NN pairwise potentials, as determined by htting 
the two nucleon scattering data, do not give the correct binding energies of systems 
with A> 3. If such potentials provide a very precise ht to the data per datum 
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~ 1), their predictions for the triton binding energy nearly coincide [6] among them¬ 
selves; nevertheless, appreciable differences continne to exist with the corresponding 
experimental valnes. The lack of binding in the three-nncleon systems with these 
models, has been attribnted to different possible canses, snch as relativistic correc¬ 
tions, non-local effects and three body interaction (TBI) terms in the Hamiltonian. 
Charge symmetry breaking (CSB) terms are also important to nnderstand the mass 
difference between and ^He, bnt give a minor contribntion to the total energy. 
Of conrse, all these effects are not independent from each other and their correct 
description is a difficnlt task. 

One of the most important aims in few-nncleon theory is testing realistic NN 
potentials, not only in the bonnd states, bnt in continnnm states too, mainly in 
those related to the A=3 systems. As a matter of fact, a series of theoretical 
investigations has been devoted to the stndy of the p-d and n-d elastic scattering 
and break-np reactions. The general statns of the field is satisfactory [7] and an 
overall agreement between theory and experiments has been reached. Among the 
possible improvements, here we shall be interested in a correct treatment of the 
Conlomb forces and in a detailed comparison of theoretical and experimental resnlts. 
To this end we have extended the method of ref.[l] to describe N-d scattering 
processes and, in this paper, onr attention will be devoted to the stndy of states 
below the denteron break-np treshold. The break-np reactions will be the object 
of a snbseqnent paper. 

We are here particnrarly interested to calcnlate scattering lengths and phase 
shifts for the varions N-d elastic scattering channels, to be compared with the 
existing experimental analysis [8]. However, snch qnantities are strictly related 
to the gronnd state energy ( for example, the donblet scattering lenght and the 
binding energy he on the so-called Philips line [9]). As already atated, the AV14 
interaction alone, as considered in ref.[l], is not adeqnate and for this reason, we 
have inclnded TBI terms in the the Hamiltonian. These terms have been proposed to 
take into acconnt all the important two pion exchange effects, in an approximate bnt 
qnantitative way. We will consider in this paper the so-called Tucson-Melbonrne 
(TM) [10] and Brazil (BR) [11] three-body potentials. 

One important aspect of the present calcnlations is the extension of the vari¬ 
ational method to scattering states. The Rayleigh-Ritz variational principle has 
a long tradition in nnclear physics and it has been largely nsed, starting from the 
first attemps to investigate the bonnd state of few-nncleon systems. On the other 
hand, a few calcnlations exist for scattering states and, among them, those nsing 
variational techniqnes, as the Kohn-Hnlten principle, are scarce. The reason is re¬ 
lated to the difficnlty of obtaining accnracies comparable with those of bonnd state 
calcnlations. 

In recent times, accnrate variational techniqnes were developed to solve the 
three- [1,3] and fonr-nncleon [12] problems, showing a fast convergence in the nnm- 
ber of channels when compared with the Faddeev techniqnes. It is natnral to ask 
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whether the situation would be similar in the case of scattering states. This is one 
more motivation for the present work. 

The paper is organized in the following way: in section 2 the construction of 
the trinucleon bound state w.f. in terms of the PHH basis functions is worked out 
and the resulting equations are solved with the inclusion of TBI terms; in the next 
section the method is extended to describe N-d scattering states below the deuteron 
break-up threshold; the basic equations which allow to determine the w.f and the 
reactance matrix 3^ are obtained in the frame of the Kohn variational principle; the 
numerical procedure used to solve those equations is shortly outlined in sections 2 
and 4 and, in the latter the results obtained for the scattering lengths and phase 
shifts are given. Finally, the merits of the approach are discussed in the last section, 
together with its possible extension to other problems. 

2. The three—nucleon bound state. 

The three-nucleon Hamiltonian is taken to have the form 

H = T + Y,Vit.j)+ WihJ.k), ( 2 . 1 ) 

i<j i<j<k 

where T is the non relativistic total kinetic energy operator and the two- and three- 
body interactions, V and IF, are explicitely included. The Coulomb p-p repulsion, 
if it is present, is contained in the two-body potential; moreover, the three particles 
are assumed to have the same mass (h^/M = 41.47 MeV fm^ throughout the paper). 
Let us introduce, in the center of mass (c.m.) frame, the following set of Jacobi 
coordinates (f, j, /c = 1, 2, 3 cyclic) 


= (rj-rfc), yi = —{vj + Yk 


2r,: 


The w.f. of the system is written as a sum of three Faddeev amplitudes. 


( 2 . 2 ) 


4' = </>(xi,yi) + </>(xi,yj) + . (2,3) 

In the above equation, each amplitude corresponds to a total angular moment um 
J Jz and total isospin TT^, therefore, if we use the L-S coupling, it can be written 
in the form 


'0(xi,yi) 

yaijk.i) 


N, 

= ^ ^a{xi,yi)yaijk,i) 

a=l 



TT^ ’ 


(2.4a) 

(2.46) 
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where yi are the moduli of the Jacobi coordinates. Each a- channel is specihed 
by the angular momenta coupled to give Aq,, and by the spin (isospin) 

(t^a) Sq (t^) of the pair j, k and the third particle i, coupled to give 
(T). The number of channels taken into account to construct the w.f. can 
be increased until convergence is reached. The antisymmetrization of the w.f. T 
requires la + to be odd; in addition la + La must be even for positive 

parity states and odd for the negative ones. 

The ground state of the three-nucleon system has positive parity and J = 1/2, 
T = 1/2 (the inclusion of the Coulomb potential gives a small T = 3/2 contribution, 
disregarded in the present work, as well as other charge symmetry breaking terms 
in the Hamiltonian). The channels allowed by such conditions are easily obtained 
and ordered of increasing angular momenta values; for example, they are 10, 18 and 
26 for £« + Tq, < 2, + Lq, < 4 and la + La < 6, respectively. 

Let us introduce, in place of the coordinates Xi^ yi, the hyperspherical coordi¬ 
nates dehned by: 


Xi = p cos(f)i, yi = psm(f)i, (2.5) 

where p is the hyperradius. The radial dependence of each a-amplitude in the w.f. 
(2.4a) is now expanded in terms of the PHH basis functions in the following way: 


$a(a:i,yi) = ^ u^{p) 


TD^Oi ? q ; 




K=Ko 


( 2 . 6 ) 


where the hyperspherical polynomials are given by [13] 

(2)p^“’''"“((/)j) = iV,^“’-^“(sin(/)i)'^“(cos(/)i)^“P^“+^/2,C-ri/2^PQg2(/)j) , (2.7) 

jg normalization factor and is a Jacobi polynomial. The grand orbital 
quantum number is given by AT = J- Lq, J- 2n, with n a non-negative integer. In 
eq.(2.6) Kq = la + La is the minimum grand orbital quantum number and Ka is 
the maximum selected value, so that the number of basis functions per channel is 


Ma = {Ka-Ko)/2 + l, 


( 2 . 8 ) 


corresponding to the maximum value of the index n plus one. When a goes to 
inhnity, the expansion basis used in eq.(2.6) is obviously complete. 

If the functions fa{xi) in eq.(2.6) are taken equal to one, the standard (uncor¬ 
related) HH expansion is recovered. Such an expansion is well suited to describe 
the structure of the system in the case of soft interparticle potentials, where a 
rather small number of basis functions is sufficient to reproduce the w.f. within 
a reasonable accuracy [14]. However, for potentials containing a strong repulsion 
at small distances, the w.f. must be accurately determined for small interparticle 
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separation values and correspondingly the rate of convergence of the HH expansion 
results to be very slow [15]. The role of the correlation function fa{x) in eq.(2.6) is 
therefore to fasten the convergence of the expansion by improving the description 
of the system when a pair of particles are close to each other. A simple procedure 
to determine the correlation functions is the one outlined in [ 1 ],where the functions 
fa (r) are taken as the solutions of the following zero-energy Schroedinger equations 


W + Vp,f3'ir) + Xp,p'{r)]fpfr) = 0. (2.9) 

/?' 

T/ 3 ^/ 3 / and Vg,/ 3 ' are the kinetic and the potential energy operators, 


T, 


P,P' 


m 


,29 ip{ip + l) 

df.2 




( 2 . 10 ) 


Vf3,(3' = < ^/3 4^ 4^ I1^/3 4^ 4^ > ’ 

and V{jk) is the NN potential; the term Xp^piir) in (2.9) is chosen of the simple 
form [1,16] 


A/3,/3'(r) = A5^exp(-7r)h/3,/3/ (2.11) 

and its role is to allow the function f/ 3 {r) to satisfy an appropriate healing condition. 
This can be achieved by taking 7 as a trial parameter, whose precise value is not 
important (the choice 7 ~ .5fm“^ is adequate [16]), and conveniently hxing the 
depth A 4 As an example, in the case of an uncoupled channel we can require that 

fp{x) = 1, when r > R, (2.12) 

where R is large with respect to the range of the potential For coupled 

channels, the condition ( 2 . 12 ) is satished by the function associated with the lower 
angular momentum value, while the other function goes to zero when r becomes 
large. 

The remaining problem is the determination of the hyperradial functions uf^{p) 
contained in eqs.(2.6); to this aim we will use the Rayleight-Ritz principle, requiring 
that the following condition to be satished 

< - E\-^ >=0, (2.13) 

where represents the change in the w.f. caused by an inhnitesimal variation of 
the functions 'u^(p). From the latter equation, it follows that 

/a(a:,) >^=0, (2.14) 
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where the subscript O indicates that the integration over the hyperangles (fti and the 
angles Xi,yi must be performed, and dVL = sin^ (picas'^ (pi d(j)idxidyi. From eq(2.14), 
after the evaluation of the spin-isospin traces and the angular integration, one 
obtains a set of second-order differential equations for the functions u'^,{p) which 
can be written in the form [1] 





dp 


'K,K 


771 

,{p) + ^EN 
n 


a,a 
K,K' 


ip) Uk'{p) 


0, (2.15) 


with a' = 1,..., Nc and K' = Kq, ..., K'^,. A numerical technique for solving the 

set of eqs.(2.15) has been outlined in ref.[l]. An important point is to reduce as much 

as possible the number N of grid points [pi,..., p^v] in the hyperradius, without 

losing accuracy in the solution. Due to the nature of the problem, the coefficients 
/ 

{X = A,B,C,N), and correspondingly the solution, strongly vary only for 
small values of p; for this reason, a grid of the type Pk+i — Pk = x{Pk — Pk-i), where 
the step-length is increased by a constant factor, results to be a convenient choice. 
As a consequence, a new variable w is introduced with equally spaced grid values 
Wk and the corresponding pk values can be obtained by the relation 

with h = p 2 — pi. For potentials without hard-core repulsion, as the ones considered 
here, we have pi = rci = 0 , and the last grid point can be chosen so that 
Pn = Pmax is about 20 -h 25 fm for bound state calculations. The set of eqs.(2.15) 
can be re-written in terms of the variable w and the corresponding generalized 
eigenvalue problem then solved by standard procedures. 

Another useful possibility is to introduce a new variable to map the inhnite 
range of the hyperradius [0, cx)] into the hnite interval [0,1]. The new variable is 
then transformed into the variable w which is used with a constant step grid. This 
is achieved by the following two transformations: 


z = 1 — exp(—i/p) 

«.= ^lnll + i(x-l)l . (2.17) 

Iny h 

where the parameters iz and x are chosen in such a way to get the most convenient 
set of Pk values. The main difference between the trasformations (2.16) and (2.17) 
lies in the fact that, when solving the set of differential equations, a boundary 
condition mut be imposed at p = Pmax in the hrst case, whereas in the second case, 
the hnal point in the rc-grid corresponds to p = cx), and for such a point the usual 
boundary condition is that the function to be determined is zero. In the calculation 
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of the A=3 bound state w.f., both the transformations (2.16) and (2.17) produce 
almost exactly the same results. However, in the study of scattering states, the 
asymptotic conditions to be satished are more easily fulhlled by the transformation 
(2.17), which will be therefore the one adopted in the next section. 

The ground state of the triton has been studied in ref.[l] by means of the PHH 
approach with the AV14 potential. The system of equations (2.15) has been solved 
by including up to 12 channels in the expansion (2.6); for the sake of completeness, 
we have extended the calculations to include all the channels with + Lq, < 4 and 
the results are displayed in table 1. By inspection of the table, one can see that the 
PHH expansion is rapidly convergent to results which are in complete agreement 
with those obtained in ref. [3], where a different technique for constructing the w.f. 
has been used. These circumstances should eliminate any doubt about a still missing 
small contribution to the binding energy, due to the use of variational bases in both 
the approaches. 

We have also studied the structure of the A=3 systems when three-body forces 
are included in the Hamiltonian of the system. Therefore, we have implemented the 
nucleon interaction by adding to the AV14 potential both the Tucson-Melbourne 
(TM) and the Brazil (BR) three-body potentials, which have been used by differ¬ 
ent groups [3,17,18] to study the three-nucleon bound state. In these models the 
numerical choice of the ttN form factor cutoff A appears to be critical, the choice 
A = 5.8/U (/U is the pion mass) checked in ref.[17], leads to an overbinding of the 
trinucleon system. Of course, it is possible to consider also other TBI effects, as for 
example those due to p-meson exchanges [19]. In order to compare with the men¬ 
tioned calculations, as a hrst choice we have adopted the same A value of ref. [17]. 
The results obtained by including up to 18 channels are shown in table 2, together 
with the estimates of the Los Alamos group corresponding to a 34 channels con- 
hguration space Faddeev (CSF) calculation [17]. The agreement between the two 
approaches is very satisfactory. Successively, the value of the cutoff parameter A 
has been changed in order to provide a triton binding energy value close to the 
experimental one. The results are listed in table 4 for the AV14 -|- TM and AV14 
-|- BR models and correspond to the values A = 5.13p and A = 4.99//, respectively. 
As mentioned in the introduction, it is necessary that the adopted potential models 
give a satisfactory description of the A=3 bound states, due to the correlation be¬ 
tween low energy scattering observables and ground state binding energy, and for 
this reason the modihed values of the A parameter previously given will be used in 
the next section. 

Finally, the ^He ground state could easily be calculated with the same TBI 
terms and with the Coulomb potential. However, it is well known that the Coulomb 
interaction alone is not able to completely explain the mass difference between ^H 
and ^He and other CSB terms must be added to the Hamiltonian [18,19]. This 
problem will no longer be discussed here. 
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3. N-d scattering below the break-up threshold 

The progress in the experimental and theoretical stndy of N-d interaction pro¬ 
cesses dnring the last decade has been noticeable. In particnlar, N-d elastic scatter¬ 
ing and break-np cross sections data are now available over a large range of energies. 
Accnrate resnlts have been also obtained for polarized nncleon and denteron tar¬ 
gets [8]. The theoretical analyses have been done nsing a variety of techniqnes 
and NN potentials [20,21], inclnding three-body forces too [20,22]. Moreover, the 
treatment of the Conlomb interaction for the pd system has been improved in a 
satisfactory way [23]. Nevertheless, the interest in this held remains noticeable, 
hrst of all becanse of the possibility of testing models and techniqnes by comparing 
the calcnlated and experimental valnes for a nnmber of important observables. For 
snch a reason, it is evident the importance of increasing the accnracy both in the 
experimental and theoretical stndies. 

The variational approach based on the nse of PHH correlated fnnctions can be 
extended to investigate scattering states and in this section the application to the 
N-d scattering below the break-np threshold is discnssed. Following the pioneer¬ 
ing work of Delves [24] for realistic NN interactions, the wave f un ction for a N-d 
scattering state will be written as 


T = Tc + • (3.1) 

The hrst term Tc mnst be snfhciently hexible to gnarantee a detailed descrip¬ 
tion of the “core” of the system, when the particles are close to each other and the 
mntnal interaction is large; Tc goes to zero when the nncleon-denteron distance 
rNd increases. As in the previons section, Tc is the snm of three Faddeev ampli- 
tndes which are in tnrn expanded in terms of the PHH basis fnnctions. The second 
term of eq.(3.1) has to describe the asymptotic conhgnrations of the system, for 
large r^d valnes, where the nnclear N-d interaction is negligible. In the asymptotic 
region the w.f. T rednces to which therefore mnst be the appropriate asymp¬ 
totic solntion of the Schroedinger eqnation. Tyi can also be decomposed in three 
Faddeev amplitndes and each one of these is written as a linear combination of the 
following fnnctions 

^Lsj(xi,yi)= 

la=0,2 

. (3.2) 

In this eqnation wi^ {xi) is the denteron wave f un ction component in the waves with 
/q = 0,2] L is the relative angnlar momentnm of the denteron and the incident 
nncleon, S is the spin obtained by conpling the spin 1 of the denteron to the spin 
1/2 of the incident nncleon. Therefore, an asymptotic state will be labelled as 
(2S-l-i)^^ and the corresponding phase shift as Slsj- 



The functions lZ\{yi) of eq.(3.1) can be taken as the regular {R) and irregular 
(/) radial solutions of the two-body (N-d) Schroedinger equation without nuclear 
interaction. The regular solution, denoted as can be written in the form 


T ( \ = 0 

{2L + l)kHCL{v) ’ 


(3.3) 


where r] = 2Me^/3h^/c and ( = kvNd are the usual Coulomb parameters, k is related 
to the center of mass energy Ec.m. = 3ti'^k‘^/4M and y = (v^/2)rjvd- The regular 
Coulomb function and the factors Cl are dehned in the standard way and the 
non-Coulomb case is obtained in the limit 0 (see also ref. [25]). The irregular 
solution, denoted as JCLiVi), has the form 


ICL{y) = (2L + i)k^+^CL{v)^^^^ (3.4) 

where Gl{j]tO ^^e irregular Coulomb function. The function IZ^iyji) of eq.(3.1) 
is taken equal to Xi(y) for X = R, whereas for A = / it does not coincide with 
JCL^y), since this would introduce a singular behaviour at rjvd = 0, which should 
be corrected by the Tc term in eq.(3.1). In order to avoid it, the function Gl(? 7 , C) 
in eq.(3.4) has been replaced by Gl(^, C)? which differs from the previous one by a 
regularizing factor. Of course, the detailed form of such a factor is not of particular 
relevance, and in conclusion the following simple form will be used 


Gl{V: 0 = (1 - exp-^-^‘‘)^+^GLiv, 0 ■ (3.5) 

The trial parameter ^ is determined by requiring that Gl tends to smoothly 
and the value ^ = 0.25fm“^ is found to be adequate. With the above dehnitions, 
the t-th Faddeev amplitude for the asymptotic wave function is written as 

f^LSj(x*, y*) = 0f5j(xi, Yi) + Yi) , (3.6) 

L'S' 


where the matrix elements Rll' give the relative weight between the regular and 
the irregular components. They are closely related to the corresponding reactance 
matrix (3ft-matrix) elements: 

= (2L + 1)(2L' + l)k^+^'+^CLCL' ^Rit . (3.7) 

By dehnition, the eigenvalues of the 3fJ-matrix are tanhj^gj. 

The internal part Tc of the w.f. (3.1) is decomposed in the three Faddeev 
amplitudes and each one is expanded in terms of the PHH basis, as it was done for 
the bound state 
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(3.8a) 


Nc 

a=l 

^a{Xi,yi) =P^“+^“/a(Xi)[ U^{p) 

K=Ko 


(3.86) 


with the conditions u^{p) 0 when p ^ oo. 

The form of the component is given by eqs.(3.2) and (3.6) also at small 
interparticle distances, therefore the internal f un ction ninst properly correct, 
snch a T/i behavionr, in that region. For this reason, the partial wave decomposi¬ 
tion of Tc (which is trnncated when doing calcnlations) mnst inclnde, hrst of all, 
those channels that are present in the asymptotic state (open channels), bnt the 
snmmation mnst be extended also to all the other important ct-channels compatible 
with the state to be described. In conclnsion, the total wave fu n ction corresponding 
to an asymptotic state will be written as 


= Yj [^^(xi, yi)flL5j(xi, Yi)] 
i=l,3 


= 5 : 


$c(Xi, Yi) + Yi) + Y Yi 


L'S' 


(3.9) 


The Hamiltonian connects states with the same parity and total angnlar mo- 
mentnm J, bnt with AL = 0,1,2. As a conseqnence, the 3fJ-matrix for J = 1/2 is 
a 2 X 2 matrix and it is a 3 x 3 matrix in all the remaining cases. So, for a given J 
two or three independent fnnctions (3.9) can be bnilt np by different combinations 
of L and S. 

The qnantities to be determined in the wave fnnctions (3.9) are the hyperradial 
fnnctions (see eq.(3.86)) and the matrix elements Rjj/ ■ To this aim we will nse the 
Kohn variational principle. This variational principle establishes that, for scattering 
states, the 3fJ-matrix elements, considered as fnnctionals of the wave fnnction, mnst 
be stationary [23] with respect to variations of all the trial parameters. Explicitely, 
these fnnctionals are given by: 


= ■'Rlt- < ^l'S'j\R\^lsj > (3.10) 

where Rf^f^i are the trial parameters of eq.(3.9). With the above dehnition of the 
operator C and nsing eqs. (3.2-4), we hx the normalizations of the asymptotic states 
as 
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< ^LSjl^l^LSJ > — < ^LSjl^l^LSJ >— 1 • ( 3 - 12 ) 

The variation of the diagonal fnnctionals with respect to the hyperradial fnnc- 
tions u^{p) is hrst performed. From the condition 

dul^Rll] =< Su^lsM^lsj >= 0, (3.13) 

and nsing the same procednre as in the case of the Rayleigh-Ritz variational prin¬ 
ciple for to the bonnd state, an inhomogeneons set of second order differential 
eqnations is obtained: 



+ (p) -J: + (p) + 


dp 


m 

¥ 




^K’ip) 


DiM ■ 


(3.14) 

The coefficients A, B, C, N have the same expression as in the bonnd state case, 
and the inhomogeneons term is given by 


RaxiP) — P 


_ rJ'OL^ Lc 


>a , (3.15) 


where, the snbscript LSJ has been omitted for simplicity. For each asymptotic state 
(2S+i)^^ two different inhomogeneons terms can be constrncted in correspondence 
to the asymptotic function with X = R or I. Correspondingly, two different 

sets of hyperradial fnnctions are obtained by solving the system of eqs.(3.14). 

Then, in order to get the optimnm choice for the matrix elements 
the diagonal fnnctionals (3.10) are varied with respect to them. This leads to the 
following set of algebraic eqnations 


E ■'« 


ss' 

LL' 


X 


S'S'‘ 

L'L" 


= Ti 


ss' 


LL' 


L”,S' 


with the coefficients X and Y dehned as 


(3.16) 




s's'^ 

L'L" 




S'L 


J + 


L'J 




L"J 


> 


5 

(3.17) 


vSS' 

^LL' 




R 

SLJ 


+ '^sljIRI^S' 


L'J 


> 


In the latter eqnations, '^lsj indicates the internal part of the wave f un ction con¬ 
strncted with one of the two solntions of eqs.(3.14), as previonsly obtained. 
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It is worth to notice that the variation of the diagonal fnnctionals gives dif¬ 
ferent conpled eqnations which mnst be satished by the off-diagonal elements of 
the reactance matrix [26]. However, the corresponding solntions are a hrst order 
variational estimates and the last term in eq.(3.10) is, in general, not zero, as it 
wonld happen in the case of the exact wave fnnction. A second order estimate for 
the 3fJ-matrix elements can be obtained by snbstitnting, in the second member of 
eqnation (3.10), the hrst order resnlts. As the reactance matrix is symmetric, the 
method shonld provide . The degree of violation of this condition 

give nsefnl information abont the accnracy of the solntion. 

4. Numerical results 

The numerical technique to solve the inhomogeneous linear equations system 
(3.14) deserves some attention. The relevant point is that the scattering state 
w.f. (3.1) must carefully describe both the regions of small and medium-large 
interparticle distances, the asymptotic behaviour having been properly taken into 
account throught the component. For this reason, it is convenient to use in 
place of p the variable w, dehned by the transformation (2.17). By using this 
new variable, and after replacing the differential operators by hnite differences, the 
system (3.14) is transformed into a set of algebraic linear equations which can be 
solved by means of standard methods. We have used N = 40 -i- 50 rc-grid points 
have been used, the last point (wat = 1 ) corresponds to p = cx) and the penultimate 
w value corresponds to pn-i ~ 65.0 fm. 

The hrst case considered is a zero energy scattering process. In this case, due to 
centrifugal barrier effects, the reactance matrix is diagonal and two physical states 
with L = 0 can be constructed: the J = 1/2 (doublet) and the J = 3/2 (quartet) 
states. The corresponding scattering lenghts are dehned as 

= - lim . (4,1) 

k^o k 

The hrst NN interaction considered is the semi-realistic Malhiet and Tjon [27] 
potential (I-III) acting only in the s-wave. For this potential the “ exact” results 
of the CSF method [28] are available, therefore a meaningful comparison can be 
done. Our results are = 0.702 fm and ‘^ap-d = 0.003 fm for the doublet state 

and '^ttn-d = 6.442 fm and ‘^Up-d = 13.96 fm for the quartet one. These values are 
extremely close to those given by the CSF method. 

Then we have considered the realistic AV14 and AV14-I-BR interactions, with 
the TBI cutoh parameter A = 5.8//. The calculated doublet and quartet scattering 
lengths are reported in table 4 as a function of the number W of channels included 
in the expansion of the w.f., together with the 34-channels CSF results of ref.[20]. 
The agreement between the two methods is quite satisfactory, and there are only 
small differences due to the contributions from higher channels not present the CSF 
approach. 
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If the value of the cutoff parameter is taken A = 4.99/U, so as to reproduce the 
correct triton binding energy, the resulting scattering lengths are again very close to 
the CSF ones given in ref. [20]. Concerning the comparison between the theoretical 
predictions and experimental data, here we only mention that for the n-d case the 
agreement is acceptable, whereas, as it is discussed in ref. [25], one must be very 
carefull, for the p-d system, when extrapolating to zero energy the experimental 
phase shifts obtained at not sufficiently low energies. 

The numerical analysis has then been extended to elastic Nd scattering states 
below the deuteron break-up treshold. For scattering energy different from zero, the 
reactance matrix, with the exception of states having J=l/2, is a symmetric 3x3 
matrix. Its six independent parameters are the three eigenvalues (tanhj^gj) and 
the three mixing parameters which are here introduced by adopting the formalism 
of Seyler [29]. As it was mentioned in section 3, the solution obtained by means of 
eqs.(3.16) gives different hrst-order estimates for the off-diagonals elements. How¬ 
ever, the second-order estimates, calculated with the help of eq.(3.10) should verify 
the symmetric condition with an improved accuracy. In order to stress this aspect 
of the problem, the hrst- and second-order estimates for the 3?-matrix elements, 
calculated with the AV14 potential, are shown in table 5 in the case J = 1/2+ and 
a neutron incident energy of 3.0 MeV, The hrst four rows of the table correspond 
to a calculation with a number of channels — S and M = 3 hyperradial func¬ 
tions per channel. For the successive four rows Nc — S and M = 6 have been used 
and the last four rows are obtained with Nc = 10 and M = 6. By inspection of 
the table, it can be seen that the second order variational calculation accurately 
provides a symmetric reactance matrix and the correction given by the mean value 
of the operator C goes to zero as the number of hyperradial functions and channels 
increase. Analogous rates of convergence have been found for all the states studied 
in the present work. 

Accurate experimantal data and systematic phase shifts analyses exist for the 
pd system [8], for incident proton energy values En=1.0, 2.0, 3.0 MeV. We have 
calculated the scattering states with the above energy values having in mind the 
two following motivations. First a detailed comparison with the experimental re¬ 
sults using realistic interactions (including also three-body terms) and an adequate 
treatment of the Coulomb repulsion, is of interest to judge the merits of the theo¬ 
retical underlying model; second, the calculated 3fJ-matrix elements may be useful 
in the case of future phase shifts analyses of accurate experimental data. 

Coming back to our calculation, the phase shifts with L > 2 can be disre¬ 
garded at the energies here considered, so only twenty independent parameters are 
contained in the cross section, i.e. thirteen phase shifts and seven mixing param¬ 
eters. The variational estimates of these twenty parameters are displayed in table 
6 at three different values of the scattering energy. For sake of completeness, the 
system nd has also been studied and the results are reported in the table. The 
contribution from the TBI terms is unimportant in all the states, with the excep- 
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tion of the J=l/2"'", a circumstance already noticed in ref. [22], The values listed 
in table 6 correspond to the AV14 potential, but when the AV14+BR predictions 
differ appreciably from the AV14 ones, they are reported in parentheses in the ta¬ 
ble. Concerning the number of channels, the value Nc = 14 for J = 1/2"'" and 
Nc = 18 for J = 3/2"'" have been used. For the other J values Nc has been hxed 
by including all open channels (giving the major contribution) and a few selected 
close channels which give a minor but appreciable contribution; in any case the 
total number of channels does never exceed Nc = 18. The number of hyperradial 
functions per channel can vary in the range M = 3 -h 6, with enough flexibility in 
the radial functions to get convergence for all the calculated matrix elements listed. 

5. Discussion and conclnsions 

The study of the scattering processes involving a few strongly interacting par¬ 
ticles is a held of large interest, but only in recent times it has been possible to 
performe accurate numerical calculations of the important related quantities. An 
emblematic situation is that one involving three nucleons and, in particular, the 
N-d interaction processes. For this system, the main progress has been realized by 
the Faddeev method, both in momentum and coordinate representations. A deli¬ 
cate aspect of such an approach is connected with the treatment of the pp Coulomb 
interaction. As it is well known, the nuclear and Coulomb potentials are expanded 
in channels with increasing angular momenta and the expansion is enlarged until a 
satisfactory convergenge for the w.f. of the system and the (part of) potential taken 
into account is found. However, the expansion in partial waves of the Coulomb in¬ 
teraction is slowly convergent and the problem of the missing contribution must be 
carefully investigated. 

On the other side, the merit of a variational treatment is that even if the w.f. is 
expandend in channels, the interaction is fully taken into account as no expansion 
of the potentials is performed. Of course, the variational tecnique must be devised 
with particular attention in order to reproduce all the relevant details of the inves¬ 
tigated structure. This might appear as easy to be satisfy, since the number of trial 
parameters can be increased and the computational facilities available at present 
allow to handle quite easily the corresponding numerical problem. However, this 
is the situation only for rather simple homework interactions and the variational 
analyses of the bound and scattering states of the three-nucleon system, with re¬ 
alistic interactions, have encountered severe difficulties. One possible solution is to 
employ adequate expansion basis for the radial dependence of each channel in the 
w.f. expansion. Since the NN interaction produces strong nucleon-nucleon corre¬ 
lations, it appears to be convenient to introduce sets of basis functions which in 
somehow take care of these correlations. A choice, investigated with success in the 
bound state of the A=3 systems, is a correlad harmonic oscillator (CHO) basis [30. 
However, in order to reproduce an exponential type behaviour characteristic of the 
large interparticle distances of the bound state w.f. of nuclear systems , a great 
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number of HO functions is necessary. With respect to this, the use of correlated 
hyperspherical harmonic basis appears more convenient, and the results obtained 
in ref.[l] for the bound state and, in the present paper, for elastic scattering states 
on three nucleons are very satisfactory. 

There are a few aspects of interest in the results presented throughout the paper 
which we wish to point out. First, there are no difficulties in treating the Coulomb 
interaction for A=3 scattering states, even above the deuteron break-up threshold. 
The inclusion of TBI terms in the nuclear interaction does not cause additional 
problems and it is not necessary to increase the number of the w.f. channels. Again 
the reason lies in the fact that no expansion in partial waves of the TBI is required. 
Concerning the numerical studies presented for the N-d elastic scattering processes, 
we have considered one of the available so-called realistic interactions, i.e. the AV14 
potential with the Brazil three-body interaction, with a value of the A parameter 
chosen to produce a correct value of the triton binding energy. However, the exten¬ 
sion of the variational method based on PHH correlated functions to other forms of 
local or non local interactions, as for example the Bonn potential, does not present 
difficulties and could be useful for precise testing of the adopted model. Indeed, the 
theoretical analyses allow for accurate evaluations of the phase shifts and mixing 
parameters in the various channels. Moreover, since there is an overall satisfactory 
agreement between theoretical and experimantal data, future improvements in the 
accuracy of the experimental data would certainly be of interest. 

We would notice hnally that the extension of the proposed variational method 
to scattering states for systems with larger number of particles, in particular A=4, 
can be easily performed but the corresponding numerical effort strongly increases. 
Interesting numerical results can, nevetheless, be again obtained. 
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Table Captions. 

Table 1. Binding energy (B), kinetic energy mean valne (T) and S'-, P and D-wave 
percentages for and ^He in terms of the nnmber of channels Nc- The 26-channels 
calcnlations of ref. [3] are reported in the last row for each case. 

Table 2. Binding energies in terms of the nnmber of channels Nc for the indicated 
potential models. The 34-channels CSF calcnlations of ref. [17] are given in the last 
row for the sake of comparison. 

Table 3. Resnlts for the triton corresponding to Nc = 18, for the two potential 
model indicated. The ttN form factor is chosen to ht the experimental binding as 
is explained in the text. 

Table 4. The donblet and qnartet scattering lengths as a f un ction of the nnmber of 
channels Nc- The 34-channel CSF calcnlations of ref. [20] are given in the last row 
for each case. 

Table 5. First order and second order estimates for the reactance matrix in three 
different approximations, as explained in the text. The central colnmn correspons 
to the evalnation of the last term of eq.(3.10) which snmmed to the hrst order 
estimate gives the second order estimate. 

Table 6. Phase shifts and mixing parameters in degrees, calcnlated for the indicated 
c.m. energies. For each energy the hrst and second colnmn correspond to the n-d 
and p-d case, respectively. The nnclear potential is the AV14 and the nnmbers in 
parenthesis correspond to the AV14+BR model. 
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Nc 

B(MeV) 

T(MeV) 

Ps'{%) 

Pd{%) 

Pp{%) 

8 

7.660 

45.551 

1.128 

8.926 

0.066 

12 

7.678 

45.645 

1.127 

8.962 

0.076 

18 

7.683 

45.671 

1.126 

8.965 

0.076 

ref. [2] 

7.684 

45.677 

1.126 

8.968 

0.076 


^He 

Nc 

B(MeV) 

T(MeV) 

Ps'{%) 

Pd{%) 

Pp{%) 

8 

7.010 

44.687 

1.318 

8.890 

0.065 

12 

7.027 

44.780 

1.315 

8.926 

0.075 

18 

7.032 

44.797 

1.314 

8.931 

0.075 

ref. [2] 

7.033 

44.812 

1.314 

8.932 

0.075 


Table 1. 


Nc 

AV14+TM 

B(MeV) 

AV14+BR 

B(MeV) 

8 

9.241 

9.153 

10 

9.306 

9.218 

12 

9.315 

9.225 

14 

9.325 

9.235 

CSF 

9.32 

9.22 


Table 2. 


Potential 

B(MeV) 

T(MeV) 

Ps'{%) 

Pd{%) 

Pp{%) 

AV14+BR 

8.481 

49.31 

0.928 

9.547 

0.139 

AV14+TM 

8.480 

49.30 

0.938 

9.255 

0.161 


Table 3. 
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Nc 

AV14 

AV14 + BR 


0>nd 

0>pd 

O'nd 

0>pd 

3 





8 

1.211 

0.980 

0.048 

-1.056 

10 

1.198 

0.957 

-0.003 

-1.139 

12 

1.196 

0.954 

-0.010 

-1.145 

CSF 

1.204 

0.965 

-0.001 

-1.136 


O'nd 

O'pd 

O^nd 

Ojpd 

3 

6.383 

13.791 

6.376 

13.830 

10 

6.382 

13.785 

6.375 

13.825 

14 

6.381 

13.781 

6.374 

13.820 

18 

6.380 

13.779 

6.373 

13.819 

CSF 

6.380 

13.764 

6.381 

13.765 


Table 4. 
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M(7V,) 

JdSS' 

^LL' 

1 — — order 

< 1 jC I ^'lsj > 

2 — — order 

3(8) 


2.778 

-0.003 

2.775 


, ^13 

0.821 

0.028 

0.849 



0.850 

-0.001 

0.849 



62.07 

3.665 

65.74 

6 (8) 


2.753 

0.002 

2.755 



0.857 

-0.008 

0.849 



0.847 

0.002 

0.849 



65.84 

0.316 

65.52 

6 (10) 


2.746 

0.001 

2.747 


, -^13 

1/2 D 2 2 

2^02 

0.854 

-0.009 

0.845 


, -^3 1 

1/2 D 2 2 
-'^20 

0.845 

0.000 

0.845 



65.88 

-0.416 

65.46 


Table 5. 
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Ee.^.(MeV) 

0.667 

1.333 

2.0 


n-d 

p-d 

n-d 

p-d 

n-d 

p-d 


-17.7 

-12.6 

-27.9 

-23.6 

-34.9 

-31.4 


(-14.6) 

(-9.54) 

(-24.0) 

(-19.8) 

(-30.6) 

(-27.3) 


-1.00 

-0.79 

-2.58 

-2.28 

-3.91 

-3.62 


(-1.00) 

(-0.78) 

(-2.57) 

(-2.28) 

(-3.90) 

(-3.61) 

Vll2 

1.04 

1.19 

1.21 

1.25 

1.26 

1.26 


(1.50) 

(1.85) 

(1.60) 

(1.73) 

(1.61) 

(1.67) 

4 c 
^3/2 

-47.2 

-37.4 

-61.3 

-53.5 

-70.5 

-63.7 

^-^3/2 

0.60 

0.45 

1.55 

1.36 

2.42 

2.20 

^Ds/2 

-1.08 

-0.84 

-2.77 

-2.46 

-4.22 

-3.91 

£ 3/2 

0.65 

0.83 

0.72 

0.79 

0.78 

0.84 

C 3/2 

-0.11 

-0.09 

-0.23 

-0.20 

-0.37 

-0.32 

%/2 

0.55 

0.53 

1.01 

0.98 

1.44 

1.39 

^D,/2 

0.57 

0.45 

1.53 

1.34 

2.38 

2.17 

'D^li 

-1.14 

-0.91 

-2.98 

-2.64 

-4.57 

-4.23 

£ 5/2 

-0.29 

-0.37 

-0.30 

-0.34 

-0.31 

-0.35 

^-^7/2 

-1.06 

-0.84 

-2.73 

-2.42 

-4.15 

-3.84 

^Pll2 

-4.54 

-3.61 

-7.66 

-6.83 

-9.51 

-8.85 

^Pl/2 

12.2 

9.30 

19.9 

17.5 

23.9 

22.0 

£ 1/2 

2.92 

2.51 

3.95 

3.48 

5.75 

4.46 


-4.51 

-3.58 

-7.57 

-6.76 

-9.34 

-8.72 

to 

14.2 

10.9 

22.6 

20.0 

26.2 

24.5 

£ 3/2 

-1.02 

-0.86 

-1.45 

-1.25 

-1.91 

-1.67 

4 ID 
-^5/2 

13.2 

10.1 

21.5 

18.9 

25.6 

23.6 


Table 6. 
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